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Abstract 

We study numerically the probability distribution of the Yang-Lee zeroes inside 
the Griffiths phase for the two dimensional site diluted Ising model and we check 
that the shape of this distribution is that predicted in previous analytical works. By 
studying the finite size scaling of the averaged smallest zero at the phase transition 
we extract, for two values of the dilution, the anomalous dimension, r/, which agrees 
very well with the previous estimated values. 



1 Introduction 



The Yang-Lee theorem provides a theoretical, and powerful, tool to study phase transitions. 
In systems without disorder (e.g. the usual 4> 4 theories or Ising models, XY model, etc) 
this theorem allows to characterize and to estimate numerically the phase transition and 
the anomalous dimension [|]. 0. 

In the disordered case (i.e. systems with a random interactions) the theorem provides a 
tool to study (and to define) the Griffiths phase (or in other words the Griffiths singularities) 
0. The Griffiths phase is a peculiar phenomenon of disordered systems. Roughly, it is a 
region above the critical temperature of the disordered system and below that of the pure 
system (for some choices of the disorder distribution this temperature could be infinite 
HI). Below the critical temperature of the pure system, which we denote T c (p = 1), but 
above the critical temperature of the disordered one, which we denote T c (p), there exist 
magnetized domains (geometrical clusters, since of course, the total magnetization is zero, 
as we are still in the paramagnetic phase of the diluted system). These domains of non- 
zero magnetization induce a complex singularity (Yang-Lee zeroes) in the free energy as a 
function of the magnetic field (Griffiths singularity ||). 

In classical statistical mechanics the Griffiths singularities are essential singularities and 
so have no effect on the static properties of the system (nothing diverges in the Griffiths 
phase, except at the critical pointQ). 

But dynamically this phase induces a slow behavior in the spin-spin autocorrelation 
functions , the dynamic of the system becomes slower than in the "usual" paramagnetic 
phase @. 

For instance, in the three dimensional spin glass case, numerically there is a change in 
the autocorrelations functions from those of the paramagnetic case (C(t) ~ t~ x exp(— at)) 
to a short range correlations (like a behavior : C(t) ~ t~ x exp(— at 13 ), (3^1) just at the 
critical point of the pure Ising model. Obviously at the critical point of the 3d spin glass 
there exists another change in the behavior of the autocorrelation function to a spin glass 
regime 0. 

In this paper we will focus our attention on the probability distribution of the smallest 
zero in the Griffiths phase and we will confront our numerical results with the analytical 
prediction of reference ||. We have obtained a clear numerical picture about the construc- 
tion of the Griffiths singularities. 

We will also extract, using the scaling of the average of the smallest zeroes at the critical 
point, the anomalous dimension of the system and we will compare this value with previous 
numerical simulations of the system ||. 



1 In the quantum case the singularities are stronger ||. 

2 It is possible to demonstrate rigorously that for Ising like models (diluted, spin glasses, etc) the 
behavior must be: C(t) ~ exp(— a(logt) d ^ d_1 ^). I thank F. Cesi for pointing this fact to me 0. 
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2 Yang- Lee singularities 



By regarding the partition function of the pure Ising model in a finite volume L d as a 
function of the variables 

P = e~ 2h , r = e^, 

where h is the magnetic field and f3 is the inverse of the temperature, Yang and Lee || [10| 
found that the complex zeroes of the partition function in the p variable lie in the unit 
circle and there are no zeroes on the real axis. Moreover in the thermodynamical limit, 
and for (3 > /3 C , the point p = 1 becomes an accumulation point giving rise to a singularity 
in the free energy. 

Near the critical point, in the paramagnetic phase, the imaginary part of the zero 
nearest to the real axis, h s , behaves 

^~(/3 c -/3) A , (1) 

and then, in the standard way, we can write down the finite-size dependence of h s at the 
critical point 

h s ~ L^ A/u \ (2) 

Using the scaling relation A = vd — (3, where d is the dimension, we can rewrite the last 
equation as 

h s ^L-^' 2 \ (3) 
Below the phase transition, in the ferromagnetic phase, the scaling law is 

h s ~ L~ d . (4) 

In the disordered case, each sample will have a smallest h, that we hereafter denote as 
h e . We will investigate numerically the functional form of the probability distribution of 
h e , that we will write as p(h e ). 

There are some analytical results about the density of the zeroes in the Griffiths phase. 
The authors of reference 0] obtain for the density of zeroes, with imaginary part hi, of a 
diluted Ising system with a proportion of spins p the following law 

/, \ ,AlogQ) 
p{hi) oc exp( — ), (5) 

as hi <C 1, which is a very weak dependence. 

It has been assumed that a cluster of size L introduces a zero, which induces the 
previous law, that scales as (see equation (f|)) 

hi = A (6) 

where A is the inverse of the site magnetization of the cluster [|J . 

It is possible to obtain a better estimate of the prefactor of in the exponential of 
the formula @ using a variational method 0] . 

The important point is that there is a finite probability to have a zero in any neighbor- 
hood of h = 0. 
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To complete this discussion we will add that at the critical point the density arrives with 
a non zero slope to the origin, in the broken phase the density at the origin is finite, and 
above of the critical temperature of the pure system the density is zero in a neighborhood 
of the origin [|J. 



3 The model and the numerical method 

The simplest disordered system is the diluted Ising model. This model describes, for 



instance, the Anderson localization [11]], and has been studied analytically (using the map- 



ping to a O(N) theory with cubic anisotropy in the limit N — > 0) [11], [12| and numerically 



1001 

The Hamiltonian of the two dimensional site diluted Ising model in a hypercubic lattice 
of size L with periodic boundary conditions is 

H e = - 2J eiej&i<Jj, (7) 

<ij> 

where < ij > denotes nearest neighbors pairs, <jj are the usual spin variables and e« 
are independent quenched noises which are 1 with probability p and with probability 
1 — p. Obviously the system will have a phase transition only if p > p c where p c is the 
percolation threshold for the d- dimensional site percolation. For instance, in two dimension 
p c = 0.592746 p| . 



There are analytical results for this model mainly by Dotsenko and Dotsenko, and 
Shalaev [[U| (DDS) using Renormalization Group techniques. There is a change in the 
functional form of the specific heat (from log \ t\ to log[l + a log \t\], where \t\ is the reduced 
critical temperature and a is a constant), but there is no change in the v exponent. This 
result must hold for a lower dilution of spins. For this weak disorder there are numerical 
results that support this picture [ |T7]j . 

But the authors of reference |§ claim that the specific heat follows the prediction of 
(DDS) but only for a lower degree of dilution, moreover they found a dependence of the v 
and 7 exponents with the dilution such that the r\ exponent is constant (we remark that 
7/z/ = 2-r?). 

The end-point of the critical line (in the plane (f3,p)), (/3 = oo,p c )Q, has critical expo- 
nents v = 4/3 and 7 = 43/18 which implies r] = 5/24 ps 0.2083 BT5 . 



The partition function for a purely imaginary magnetic field, ih, in a d- dimensional 
lattice of size L is 

Z{{3, h)=J2 exp(/3 E Wi + ih E (8) 

[a] <ij> i 

By defining M = J^i^i, the total magnetization of the system, we obtain 

Z{p, h) = (cos(JiAf)) + i{sm(hM)), (9) 

where the average ((••)) is taken with Z(/3, h = 0), i.e. a real measure. In the paramagnetic 
phase all the odd moments of the magnetization vanish, which implies (sin(hM)) = and 
the only singularities of the free energy (logZ(/3, h)) will arise from the zeroes of {cos(hM)) . 

3 This is the two dimensional site percolation phase transition. 
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This is the scenario for the pure systems. In the diluted case we need to replace cr, by 
ejCTj and so each samples will have its own smallest zero (h e ). The averaged values over all 
the samples, h e , should follow the previous finite-size scaling relation (§) at criticality. 

4 Probability distribution of the Yang-Lee zeroes in 
the Griffiths phase 

To check the analytical form of the probability distribution, p(h e ), of the smallest Yang-Lee 
zeroesf] we have done numerical simulations with (3 = 0.52 and p = 0.889 which is inside 
of the Griffiths phase. We used the Wolff algorithm |T^[ and we simulated the sizes L = 4 
(8000 samples), L = 8 (15000 samples), L = 12 (2200 samples) and L = 16 (3926 samples). 
The results are shown in Figure 1. 




Figure 1: Probability distribution of the smallest zeroes for (right to left) L = 8, L = 12 
and L = 16. 

We expect that the minimum value of h e , at fixed L, should be due to the sample with 
all the points filled (i.e. a pure Ising model of size L). We find that the minimum smallest 
zero, as a function of the lattice size, follow the rule 

Lmin 1 " 5 ( 1 ) (IQ) 
n t £ L93 ( 4 )' ^ ' 

4 Obviously in p(h) are all the possible zeroes, but as we are interested in the h <§; 1 regime then 
p(h)^p(h). ' 

5 We remark that for this dilution the phase transition is at (3 — 0.5380(3) and the phase transition of 
the pure model is at f3 = \ log(l + V2) = 0.44069. 
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where we have fitted using 4 < L < 12 with x 2 /DF = 0.2/1 (DF means degrees of freedom). 

For L > 16 lattices the previous fit ( [370] ) does not hold. This discrepancy comes from 
the fact that the number of samples that we need to pick up this minimum value is larger 
than the number that we have simulated. 

Simulating directly the pure Ising model we find that the smallest zero (simulating up 
to L = 32), as a function of the sizes, at (3 = 0.52 (ordered phase of the pure model) 
behaves 

^in = 1^5(4) 

"'pure £1.97(1) ' ^ ' 

following very well the law (HI). The agreement with eq. ( |I~0| ) is also very good. 

We have also fitted the mean value of the probability distribution as a function of L 
(for L=8, 12, 16, 64, 128, 192) in the diluted case and we have found that the numerical 
data behave 

K = 0.0041(1) + 1.67(2)L _L84(1) , 

with a very good x 2 /DF = 4.9/4. The samples used for L > 64 are written in Table 1. 

We plot in Figure 2 the head (i.e. the region of lower values of h) of the probability 
distribution for the L = 8 lattice in the variables (1/ h t , log p(h e )) in order to check the 
formula (||). 

We see two different regions that we mark with two linear fits. The first region (left part 
of the figure) has a slope (— 0.11(3))[] which agrees, is a two standard deviation, with the 
naive theoretical prediction ((logp)A = (1.5(1) x log|) = —0.18(1)), where we have used 
for A the numerator of the fit fllPf) . The second region decays with a behavior compatible 
with the equation ([5]) but the slope is wrong (slope= -0.94(4) )Q. We think that this decay 
is due to a finite-size effect (the lattice size is 8) and hides the decay with the "naive" 
slope (w -0.18). 

Bray |EJ shows that the real slope (in absolute value) has as upper bound the "naive 
value" (0.18). Our numerical results go in this direction. In particular as T — > T c (p) + the 
real slope, in absolute value, goes to zero however the "naive" value will clearly be different 
from zero. 

Hence, the numerical picture is as follows (we remark that we are in the Griffiths 
phase): we have a narrow probability distribution with its mean value having a non zero 
thermodynamic limit. But the minimum value of this probability distribution follows the 
law of the pure Ising model in the ferromagnetic phase so that goes to zero and introduces a 
singularity in the free energy. We have seen this behavior when simulating a large number 
of samples up to L = 12. Using a very large number of samples could be possible to 
continue this result to large lattices (L > 16). 

We will see in the next section how, at the critical point, the mean value of the smallest 
zeroes goes to zero following a power law. 

6 Result of a least square fit using the points second to fourth in figure 2 (seen left to right). 

7 Using the points sixth to ninth in figure 2 (seen left to right). 
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Figure 2: Head of the logarithm of the probability distribution (modulo a normalization 
factor) of the smallest zeroes for L = 8 as a function of l/h e . 



L 


p = 0.889 


p = 0.75 


64 


100 


100 


128 


40 


40 


196 


30 


40 


256 


30 


40 



Table 1: Number of samples simulated for different sizes and dilutions used in the numer- 
ical simulations of sections 4 and 5. 

5 Scaling of the Yang-Lee zeroes at T c 

At T c we have performed numerical simulations using the Wolff algorithm with two degree 
of dilution, p = 0.889 and p = 0.75, and lattice sizes L = 64, 128, 192 and 256. We report 
in Table 1 the number of samples used. 

We have used the values of the inverse critical temperatures (f3 c (p)) reported in reference 
|| 0. We will also compare the results of this reference with our result for the rj exponent. 

We have measured the susceptibility, 

where (■■) is the average on the disorder and ((■■)) is the thermal average. We have also 
8 i.e. p c (p = 0.889) = 0.5380(3) and /3 c (p = 0.75) = 0.772(1). 
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V 




V 


l/v 


A/z/ 


V 


0.889 


1.72(1) 


0.279(14) 


1.75(2) 


1.873(13) 


0.254(26) 


0.75 


1.72(3) 


0.28(3) 


1.76(3) 


1.89(2) 


0.22(4) 



Table 2: Results for the critical exponents. The first column is the proportion of spins. 
The next two columns are the critical exponents reported in reference [8]. The fourth and 
fifth ones are our estimates of 7/z/ (as control, calculated as Xmax ~ L 1 ^) and A/z/ (using 
the scaling of the zeroes) respectively. In the last column we report r\ calculated using the 
scaling relation 77 = 4 — 2 A/z/. 

measured (cos(hM)) in every sample to calculate the zeroes. 

We obtain h e , the smallest zero for each sample, and then we calcule the mean value, 
h e . The error is estimated using sample to sample fluctuations. We plot the finite-size 
scaling in Figure 3 and Figure 4, for p = 0.889 and p = 0.75 respectively, with our best 
power fits drawn as a line (fifth column of Table 2). We report the numerical values of the 
fit (also for the susceptibility) in Table 2. The second and third columns of the Table 2, 
are the estimates of reference || for 7/z/ and 77, obtained as 2 — 7/z/. 



v 



-,■[: 
O 



log L 



Figure 3: Mean value of the smallest Yang-Lee zeroes, < h e >= h e , against the size, in a 
double logarithmic scale for a dilution p = 0.889. The straight line is the power law fit 
reported in Table 2 and text. This fit has a % 2 /DF = 0.11. 

Table 2 shows that our values of 7/z/ are in the errors with those of reference |§ (we 
perform this as check) and this also holds with our estimate of rj using scaling of zeroes. 
The results are compatibles with rj = 0.25 on the critical line. 



8 



4.5 5 5.5 

log L 



Figure 4: Mean value of the smallest Yang-Lee zeroes against the size in a double logarith- 
mic scale for a dilution p = 0.75. The straight line is the power law fit reported in Table 2 
and text. This fit has a xV DF = °- 28 - 
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6 Conclusions 



We have investigated the Griffiths phase by studying the behavior of the probability dis- 
tribution of the smallest Yang-Lee zeroes. We have obtained a clear numerical picture of 
the finite-size construction of these singularities. We have also confronted our numerical 
data with previous analytical results (4| and the agreement is very good. 

In the second part of this paper we have shown that the study of the smallest zeroes is 
very useful to estimate accurately the anomalous dimension of the system. 

We have extracted one critical exponent of the system, rj, which agrees with the ana- 
lytical predictions and with the numerical results. We need to calculate the second one in 
order to fix the universality class of the Hamiltonian. A possible calculation, in the line to 
seek complex singularities, is the study of the Fisher zeroes (TE| . This study will point out 
the thermal critical exponent v [|2(J and clarify if it depends on the proportion of spins or 
not. 
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